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Abstract 



The anomalous scaling exponents Cn of the longitudinal structure functions 
Sn for homogeneous isotropic turbulence are derived from the Navier-Stokes 
equations by using field theoretic methods to develop a low energy approxi- 
mation in which the Kolmogorov theory is shown to act effectively as a mean 
field theory. The corrections to the Kolmogorov exponents are expressed in 
terms of the anomalous dimensions of the composite operators which occur 
in the definition of Sn- These are calculated from the anomalous scaling of 
the appropriate class of nonlinear Green's function, using an uv fixed point 
of the renormalisation group, which thereby establishes the connection with 
the dynamics of the turbulence. The main result is an algebraic expression 
for (^n, which contains no adjustable constants. It is valid at orders n below 
g~^, where 5* is the fixed point coupling constant. This expression is used to 
calculate for orders in the range n = 2 to 10, and the results are shown to 
be in good agreement with experimental data, key examples being ^2 = 0.7, 
C3 = 1 and Ce = 1-8. 

I. INTRODUCTION 

The study of homogeneous isotropic turbulence has as its aim the derivation of the 
statistical features of small scale velocity fluctuations at high Reynolds numbers, based 
on the assumption that they exhibit universal characteristics independent of the form of 
the large scale flow structures [1-3]. A key quantity of interest is the longitudinal velocity 
increment, v+ — where v± = Vi{x ± r /2,y, z,t), the velocity component Vi and the 
separation distance r both being in the same direction, here the x-axis. An empirical fact 
is that its nth order moment, the longitudinal structure function <S'„(r), defined by 

Sn{r) = (K - V-T) , (1) 
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exhibits multiscaling. That is, the exponent Cn, defined by the scahng relation 

Sn{r) ~ (2) 

is a nonhnear function of the order n. This behaviour is not explained by the classical 
turbulence theory of Kolmogorov [4] which yields a linear dependence 

= I (3) 

Moreover, the amount by which (n differs from called the anomaly, has proved stub- 
bornly resistant to attempts at quantitative explanation [1-3,5,6]. The obstacle to progress 
with the theory is the strong nonlinearity of the governing Navier-Stokes (NS) equations. In 
this paper, our aim is to show how modern statistical field theory can be used to overcome 
this difficulty and provide theoretical predictions for (^n, which agree well with turbulent fiow 
data. 

The idea that statistical field theory can be brought to bear on the problem of turbu- 
lence is not itself new. Indeed, interest in describing turbulence in terms of the underlying 
functional probability distribution of the velocity field, together with its corresponding gen- 
erating functional W, has a long history [5,6]. But such work has suffered from the weakness 
of relying on conventional perturbation theory to effect closure of the statistical hierarchy, 
whereas it is widely believed that a non-perturbative treatment is necessary, because the 
NS equations lack a small parameter. Consequently, progress with this approach has been 
disappointing. 

The question is whether we can find a middle course, which avoids the limitations of 
conventional perturbation theory, while not demanding an intractable non-perturbative ap- 
proach. Here we explore the possibility of formulating a more efficient perturbation theory 
by developing a zero-order solution which already accounts for the dominant nonlinear inter- 
actions, in an attempt, as it were, to deplete the effect of the nonlinearity. We shall do this 
by adopting a more general quadratic form in W in place of the viscous form which arises 
naturally. The modified quadratic form is determined self-consistcntly from the NS nonlin- 
earity using the linear response function and the energy equation. In the inertial range, it 



leads directly to the Kolmogorov distribution, after allowing for the kinematic effect of the 
sweeping of the smaller scales by the larger ones. The difference between these quadratic 
forms then appears as a perturbation, which, as we shall see, is not critical, provided that 
the force spectrum function is non-zero only at small wavenumbers and yields a finite input 
power. 

Having incorporated the dominant nonlinearities which are responsible for the turbulence 
energy cascade into the zero-order solution, what is then lacking is the effect of the fluctuating 
dissipation rate, which is the well-known defect of the Kolmogorov theory [5]. In this 
approach, the perturbation theory is then, in effect, only required to accommodate the 
residual coupling associated with these fluctuations, which are directly responsible for the 
anomalies. The fact that the anomalies are small, and associated with a weak residual 
coupling, provides good reason to expect that a small expansion parameter might emerge, 
thereby rendering the problem accessible to perturbation theory, essentially by means of a 
standard loop-expansion of the generating functional. 

Although the use of the modified quadratic form as an initial approximation would appear 
to be an attractive option, providing a sound physical basis for the approximate evaluation 
of the generating functional, it poses severe technical problems, the most significant being 
the occurrence of divergences at higher orders in perturbation theory, due to the incomplete 
representation of the large scale flow. These divergences are of two types: power divergences 
(including power x logarithmic), which are associated with the sweeping and pure loga- 
rithmic divergences, which describe the cascade process. On the other hand, statistical field 
theory [7] , provides the mathematical techniques needed to compensate for such divergences, 
in the form of the well-known processes of resummation and renormalisation. In particular, 
renormalisation [8] provides a procedure whereby the scale invariance in (Q) can be recovered 
from a divergent theory, yielding the exponents in terms of the anomalous dimensions of the 
composite operators appearing in (|I]), which we can calculate from the appropriate nonlinear 
Green's functions. The modified quadratic form itself follows uniquely from the requirement 
for the absence of non-renormalisable terms, after renormalising the basic parameters of W 
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and allowing for sweeping. 

Fortunately, sweeping effects do not pose an insuperable obstacle, notwithstanding that 
the initial formulation is Eulerian. Indeed, we show that the power divergences associated 
with sweeping can be removed by introducing a single sweeping interaction term, which 
can be derived from the generating functional itself using a random Galilean transforma- 
tion of the velocity field, having an rms convection speed which is calculated from the NS 
nonlinearity. The application of this transformation does not, of course, affect the values 
of S'„(r) and, thus, enables the straining interactions which determine the spectrum to be 
separated from the background of sweeping convection, yielding, in effect, quasi-Lagrangian 
forms. In this way, as we shall show below, it proves possible to demonstrate multiscaling 
and calculate the anomalies of the structure functions accurately. 



II. THEORETICAL FOUNDATIONS 

Our starting point is the NS equations, describing flow in an incompressible fluid of unit 
density, velocity v, kinematic viscosity V and pressure p, and driven by a random solenoidal 
stirring force f , which are 

_ + vVv = -Vp + z/VV + f, (4) 
at 

and 

divv = 0. (5) 

Suppose that 

v(x) = V(x|f) (6) 

is the solution of (^) at the space-time point x = (x, t), corresponding to a force f (x), which 
has a Gaussian probability distribution V{f). Then the generating functional W for the 
correlation functions of the velocity field can be written as the functional integral 
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W = J exp{S)V{f)Vf, (7) 

where the source term is given by 

S{3) = J 3{x) ■V{x\i)dx, (8) 

and the correlators follow by functional differentiation with respect to the source field J(x). 
Given that we cannot obtain an explicit expression for the solution the crux of the 
problem is how to approximate with the accuracy required to calculate the (n- In our 
approach, as indicated above, we prove that the Kolmogorov theory can be used effectively 
as a mean field theory in a saddle-point evaluation of (0), and that this leads to an expansion 
which has a genuinely small coupling constant. 

Within the context of a field theoretic interpretation of (^, each term of the binomial 
expansion of (|l]) must be regarded as an operator product of the usual Wilson type, (see eg 
[7,8] ). Correspondingly, the powers of v± must be treated as composite operators, which, 
in accordance with standard procedures [7], must be generated from W by independent 
sources. Here, our aim is to limit the composite operators that need to be allowed for to 
those which appear explicitly in the definition of Sn{r), as given in ([1|). To this end, we 
define a set of longitudinal composite operators Os{x), for s = 2,3, 4, by 

Os{x)=v,{xy/s\, (9) 

which we generate from W by adding to the source term (P), the additional term 

- ^ / tsix)Osix)dx. (10) 

We also need to include in the definition of W a means of establishing the vital link 
between the time- independent definition of S'„(r) and the dynamics of the turbulence. This 
requires the introduction of a dynamic response operator, which we define to be the func- 
tional differentiation operator 



Its inclusion in the definition of W adds a final source term to S, given by 

Ja{x)J-'a{x)dx, (12) 
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where summation over repeated vector indices is implied here and below. 

The terms (|^), (|IU]) and (|T2|) together constitute the full source term for (|^) which 
becomes, therefore, 

S{J,J,ts) = / {Jaix)Vaix) + Jc{x)J^a{x) - ^^^t s{x)0 s{x)}dx, (13) 

s 

and this completes the definition of W. Thus, (0) and ([13D provide the foundation of our 
approach to the calculation of Cn- However, before we proceed with this calculation, we need 
to cast W into a conventional field theory form, and introduce the modified quadratic form. 

A straightforward method of transforming (|^) into a conventional field theory form is to 
replace 'P(f ) by its functional Fourier transform and then integrate over f . This is the stage 
at which we make explicit use of the NS equations. Essentially, to effect the transformation, 
we change our perspective by replacing the velocity field V(f) generated by the force f,by 
the force F(v) needed to excite a particular realisation v of the flow field. The operator ( pH]) 
is then replaced by an equivalent conjugate vector field v. 

To carry out this transformation, we work in the Fourier domain setting 

v(x) = j exp^ik ■ x)\^(k)Dk, 

where k denotes (k,^^;), so that ■ x = cjt — k ■ x, while Dk = dudk/ {2tt)^. Then, from (^ 
and d^), we have 

F^k, v) = GoCky'vaCk) - ^(27r)^P„^^(k) j vp{p)v^{q)5{p + q- k)DpDq. (14) 
The notation here is the following. Go{k) is the zero-order approximation to the response 



function G{k) defined below in (|63) and (B^), namely 



where 

Pai3-y{k) is the NS vertex defined by 

(k) — kpPa^{k.) + k^Pai3{k.). 

where 

Next we write the Gaussian distribution of f in the form 



P(f) = ATexp |-1 y U-k)h{kY^P^p{\^)fp{k)Dk^ , (16) 



for which the corresponding force covariance is 

'faik)m) = {27T)'6Ck + i)h{k)Pap{k), 



where the force spectrum function h{k) is an arbitrary function which is assumed to be 
peaked near the origin so that the power input / h{k)d\^ is finite. We now change the 
functional integration over f in (|^ to an integration over v by means of the transformation 
v(A;) = V(A;|f), and substitute the representation 

p(f) =M j exp|-i j v^{-k)h{k)P^p{}i)vp{k)Dk + i j v^{-k)f^{k)Dk^Vw, (17) 

Since the Jacobian only contributes an unimportant constant, we get 

iy(J,J,t,) = y exp -L(v,v) + 5(J, J,i.) ^^vDv, (18) 



where 



L(v, ^) = \j v^{-k)h{k)P^p{\^)v^{k)Dk - I j v^i-k)FXk, ^)Dk, (19) 



while the source term flTSf) becomes 



S{3,3,Q = I {U-kMk) + U-k)va{k) -Y.^s{-k)0,{k)}Dk, (20) 

^ s 

The expression ( [TsD casts W into the form of an MSR type functional integral [9] . 

Now the quadratic form appearing in (|19D does not provide a good initial approximation 
for inertial range scaling because, of course, it merely describes the viscous decay of an 
externally driven random flow, with no account taken of the nonlinear interactions. It is 
thus essential in developing an expansion theorem for (|18[) to introduce a more appropriate 
quadratic form. Now the general theory of quadratic forms in a Hilbert space indicates that 
we can introduce at most two functions. These can be taken as an apparent force spectrum 
Do{k) and an effective micro timescale To{k), which are related to the energy in wavemode 
k, Q{k), by 

Q{k) = ro{k)Do{k). (21) 

The modified quadratic form in L(v,v) is then obtained, firstly, by replacing h{k) with 
Doi^k) and, secondly, by replacing the viscous timescale T^{k) by the the effective timescale 
To{k), so that the viscous propagator (|I5|) in (|T^) is replaced by 

ICU + To{k) ^ 

Thus, we now have in place of ( p!9[ ) 

L(v,v) = ^J d^i-k)Doik)P^f,{k)vfsCk)Dk-t J v^{-k)F^{k,w)Dk (22) 

in which D^lk) and TQ{k) are, as yet, unknown functions to be determined in an appropriate 
way from the energy equation and the linear response function. The idea that one should 
replace the viscous quadratic form by a modified form was suggested originally in [10], where 
it was used in conjunction with a variational principle based on an entropy functional,but 
recent work [11] has shown that this approach contains an arbitrary element. However, 
we shall not need to invoke any additional principle, because we shall be able to deduce 
the modified quadratic form in a self-consistent way from the 1-loop expansion, as we have 
already indicated. 
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The introduction of the modified quadratic form as a basis for an expansion theorem for 
(p!8|) requires the inclusion of the difference terms as perturbations, which contributes an 
additional term to L given by 

^L^ = \j U-k){h{k)-Doik)}P^p{k)dp{k)Dk-^ j v^{-k)Mk)-' -To{k)''}vXk)Dk. 

(23) 

These terms have the same form as the counterterms introduced below in (pTf ) to accommo- 
date the pure logarithmic divergences but their role, as we shall see, is not critical as regards 
calculating the inertial range exponents. 

The derivation of the functions -Do(fc) and To(fc) occurring in the modified quadratic 
form entails a detailed discussion of sweeping convection, the structure of the Feynman 
diagrams associated with the loop expansion of W and the establishment of the condition 
for the absence from the linear response function of non-renormalisable terms. We shall defer 
detailed discussion of these topics until Sections |V11| and |V111| and, meanwhile, proceed with 



the calculation of the anomalous exponents by anticipating their forms, which, in the inertial 
range, are 

Do{k) = Dok-\ (24) 

and 

To{k)-' = (25) 

Clearly, these forms imply that the zero order solution behaves in the inertial range as if the 
fiuid were stirred by a random force with a k~^ force correlation spectrum and responds to it 
with a Lagrangian time scale oc /c"^/'^. Thus, they lead to the Kolmogorov distribution. We 
shall explain how this result follows from the generating functional in Section |V111| . The 
advantage of this approximation is that it achieves a prime requirement of any efficient 
perturbation theory, which is a zero-order approximation that already closely approximates 
the desired solution. 
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On the other hand, as we have indicated, the resulting perturbation theory yields diver- 
gences at higher orders. But these divergences can be handled by standard renormalisation 
procedures. Fortunately, as regards the calculation of C„, we need consider only logarith- 
mic divergences. As discussed above, this is because the power divergences represent the 
kinematic effect of the sweeping of small scales by larger scales. Indeed, as we shall show 



in Section [VI I|, such terms are precisely those which can be generated by applying a ran- 
dom Galilean transformation of the velocity field to W. Consequently, they can be cancelled 
by introducing the appropriate vertex into W, yielding quasi-Lagrangian approximations. 
Hence, from a purely practical calculational point of view, the effect of sweeping can be 
removed from the calculation of (n simply by discarding power divergences. We are then 
left with the logarithmic divergences, which we can sum by renormalisation group methods. 

Thus, an important implication of using the modified quadratic form as an initial approx- 
imation for the calculation of (n is that renormalisation becomes a necessary preliminary. 
So we need to identify the counterterms which arise in W under renormalisation and obtain 
the transformation rule which connects the bare and renormalised generating functionals. 
Renormahsation is applied to the viscosity and force constants appearing in and 
in the usual way by introducing renormalisation constants Z,^ and Z^, which relate their 
bare values z/q and Dq to their renormalised replacements u and D by 

= vZy and Dq = DZd- (26) 

This generates counterterms in ([T9D for the elementary fields (v and v) given by 

ALef = -AZ^i J v^{-k)T{k)-^Vo,{k)Dk (27) 

+AZd\ j Vai-k)D{k)P^p{k)daik)Dk, 

where we have defined renormalisation constant increments by 

AZy^Y) = Zy^D — 1. 

The additional renormalisation which must be applied to the composite operators also 
takes the standard form 
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{Os)b = Z,{Os)r. (28) 
The corresponding counterterm is obtained by substituting ( pBD in pO| ) to get 



where 

AZs = - 1. 



We conclude this section by giving the transformation which relates the generating func- 
tional of the bare correlation functions Wb to its corresponding renormalised form Wr. 
To provide a convenient means of handling the dependence of the correlation and response 
functions on the dimensional parameters z/q and Dq, we rescale V and f by introducing bare 
fields defined by 

1 



and 



with bare frequency 



f(fc) = (:^yf^(k,^s), 



UJb = —■ 

^0 



These bare fields preserve the form of the NS equations, apart from explicitly introducing 
the non-dimensional coupling constant, defined by 



90 = (29) 



in which the appropriateness of the numerical factor will appear later from the loop- 
expansion of W. 

Under the renormalisation (^), the bare fields are replaced by renormalised fields, to 
which they are related by 
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z 



and 



VB(k,cUB) = ( ^ ) ^r{^.^r) 



fB(k,u;B) = (^^^ f/?(k,u;H), 



where 



UJ 

= —■ 

V 



These relations follow from two requirements. First, the form of the NS equations (0) must 
again be preserved, with the new constants v and D resulting in a renormalised coupling 
constant 



Second, we have to satisfy the crucial requirement that 'P(f), as given in (|T6D, remains 
invariant under renormalisation. Indeed, satisfaction of these conditions implies the desired 
relation between Wb and Wr, which from (|T8|) and pO), is readily found to be 



1 / 7 \ ^/"^ 1 / 7 \ 7 / 7 \ 



The foregoing provides the basis of our calculation of C„, which involves the following 
four stages. First, we use (|31|) and the binomial expansion of (|l]) to develop a short distance 
expansion for S'„(r), by substituting an operator product expansion (OPE) [7,8] for each 
term, based on the operators @. As shown in Section |nT|, this yields the scaling of Snir) in 
terms of uv fixed point values of standard RG functions. The second stage of the calculation 
is to demonstrate that the required uv fixed point of the RG actually exists, and then to 
deduce the corresponding fixed point coupling constant g^. This is done in Section |rV| by 
considering the renormalisation of the linear response function, using the renormalised func- 
tional in the form obtained from (113). The third stage is to calculate the specific fixed point 



RG parameters which give the anomalous component of C„. To do this, we have to consider 
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the renormalisation of appropriate nonlinear Green's functions involving the composite op- 
erators defined in (||). These are identified and evaluated in Section |V|. Having calculated 
the anomalous scaling exponent — r„p of the pth term in the binomial expansion of Sn{r), in 
the fourth and final stage of the calculation, we derive a simple algebraic expression for (n 
by maximising r„p, with respect to integer values of p, and subtracting this maximum from 
the Kolmogorov value (|^). The results obtained for (n are presented in Section where 



they are shown to be in good agreement with experimental measurements at all orders for 
which reliable data exists. Finally, the mathematical proofs, deferred during the calculation 



of the exponents, are presented in Sections |VI1| - |D^ , and comprise: (a) the demonstration 



that sweeping effects can be eliminated by means of a random Galilean transformation of the 
velocity field; (b) the derivation of the modified quadratic form from the 1-loop expansion; 
and (c) the derivation of the dominant terms of the OPEs. 

III. THE STRUCTURE FUNCTION EXPANSIONS 

In applying the OPE technique to (j^) the first point to appreciate is that the orders 
n = 2, n = 3 and n > 4 require separate treatment. The factor distinguishing 5*2 and 
from the higher order Sn is that the latter involve composite operator products, whereas 5*2 
and 5*3 do not. Also, 5*3 is exceptional in representing a transition at which corrections to 
the Kolmogorov exponents change from positive at n = 2 to negative at n ^ 4, with 
no correction occurring at = 3 in accordance with the known exact scaling law, which is 
verified, within the present framework, in Section |VIII| . This sign change is caused precisely 



because composite operator products appear in Sn when n ^ 4. 

We begin, therefore, with the relatively straightforward case of 5*2. According to (0), we 
have 

S,{r) = 2 ((v^) - {v^v_)) , (32) 
which shows that the scaling of 5*2 is determined by the behaviour of the operator product 



as r — * 0. The form of its OPE is established in Section after the necessary 
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mathematical apparatus has been set up. Its proof is given there to the accuracy of the 
calculation, ie up to and including terms of order g"^. We shall show that the operators which 
appear in its OPE are: (a) the unit operator /, with constant coefficient -E/S, where E is the 
turbulence energy; (b) the dominant longitudinal quadratic composite operator 02(x), which 
gives the leading scaling behaviour; and (c) subdominant operators including all transverse 
operators and the longitudinal higher order composite operators Os{x). However, we shall 
only be concerned with the dominant operators and so we write the expansion as 

v+V- = ^EI + C2{r)02{x) + ..., (33) 

where the dots indicate the additional subdominant terms. The scaling behaviour of this 
operator product can be found in the usual way from the RG equation satisfied by the 
leading Wilson coefficient C2{t) [12]. 

We start by considering an arbitrary equal time correlation function of order Z, given by 

Ha^,„ai{xi,...,Xl) = {Vc,-,{Xi)...VaXxi)) ■ (34) 

If we insert ( ^31) into this correlation function, we get 

Hai...ai{xi, ...,xi,x + -i,x- -?) = — iJ^, (^^1 , + C2{r)Q''^^ {xi, ...,xi,x) + 

Z Z 6 

(35) 

Is) 

where, in general, Qai...ai is the inserted correlation function defined by 

Qal..ai {xi,...,Xux) = {Va,{xi) . ..Vai{xi)Os{x)) , (36) 

and f is a unit vector along the x-axis. 

We can deduce the RG equation satisfied by the Wilson coefficient C2 in ( p3D from (p5D, 

(2) 

given the RG equations satisfied by Ha^^ a^ and Qai...ai- To obtain the latter, we need the 
transformation rule for the equal time generator of these correlation functions, which we 
denote by Vr(")(J,t,).This follows in a straightforward manner by taking time independent 
sources in (|3T|) , and integrating with respect to ujb and ujr, with the J dependence, which 
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is irrelevant here, suppressed. To simplify the result, we shall anticipate the fact, which we 
demonstrate in Section fV|, that 

Zd = Z,. (37) 

We then get 

W;^\3,ts) = wj^\j,Zsts). (38) 

According to this relation, the bare and renormalised forms of H^^ ^ are equal. Hence, 
when we change the renormalisation scale, which we denote by /x, the Fourier transform of 
Hai...ai changes according to the RG equation 

= 0, (39) 
where V is the standard RG operator defined by 



with 



Pig) = /i^- (41) 



d/j, 

In the case of Qa}...ai, we obtain from (^) and (|38|) the relation 

iQS...jR = ZsiQS...jB, 

which leads to the RG equation 

= 7.QS....„ (42) 
where 7<j is the anomalous dimension of Os given by 

7s = /x^logZ,. (43) 



For ease of notation, we have dropped the suffix R in the RG equations ( P^ ) and (^), since 
we shall always be dealing with relations between renormalised functions. 
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We now apply the RG operator (1^) to the Fourier transform of (|35|), and make use of 
]) and (pp, to get 



= {VC2 + 72C2)Q^S...a, + --- ■ (44) 

(2) 

As this equation holds for arbitrary Qa/.-.a;, it follows that 

VC2 = -72C2. (45) 

which is the RG equation satisfied by the Wilson coefficients in (pSl). 

The standard solution of this equation, corresponding to an uv fixed point [12], now gives 
for the leading term of (|33D the scaling behaviour 

C2(r) ~ r2/3-72, (46) 

where the star denotes the fixed point value of (^5]). This result, in conjunction with ( P^j ) 
and P3|), yields the scaling exponent for S'2(r), namely 

(2 = 1 + A2, (47) 

where 

A2 = -72- (48) 

We shall calculate A2 in Section [V]. 

Consider now the general case for even orders n = 2m > 2. Introducing the general 
composite operator product 

Ass'ix, r) =Os[x + Os' - , (49) 

and taking advantage of the isotropic symmetry, we can write the binomial expansion of (|l]) 

as 

/ m-l \ 

{r)=n\/2j2 i-r (x, r) + (-)'" A„,„ (x, r) \ . (50) 

\ P=0 / 
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We can identify the dominant term of the OPE of A„_p^p by factoring out the product 
and using the fact that, by (|33|), its expansion begins with the unit operator. We 
will justify this process in Section This implies that the OPE of A„_p_p itself takes the 
form 

An-p,p(£, r) = Cp^rn-pir)02(m-p)ix) + (51) 

where again the dots indicate subdominant terms. Substituting (^) in (|50D , we get 



m— 1 



S^r) = n\l2 ^(-)PCp„(r) (02(™_p)(x)> + (-)"^C„,„(r) I + . . . , (52) 
I P=0 J 

the averages of the composite operators being independent of x for homogeneous isotropic 

turbulence. 

To find Cn from this expansion, we have to determine which term or terms on the right 
hand side yield the negative correction of maximum magnitude to Cjf°'- As before, this is 
deduced from the RG equation for the Wilson coefficient Cp,s, which we derive next. 

We begin by inserting (^) into the general correlation function ([3^ ) to obtain the general 
inserted correlation function 

According to (^) its bare and renormalised forms are connected by the relation 

(ss') 

from which it follows that Rai.'.ai satisfies the RG equation 

1^Ri% = {ls + ls')Rt%. (54) 

Next, we insert the expansion ([5l| ) into the general correlation function (|3^) , and use the 
definitions (|36|) and (0), to get 

R-ai-.^a!^ ^^1) ■■■iXl^X + — X — — = Cp^m-p (''") Qail.ai (-^l) ■■■^ ^) ■ 



We then apply the RG operator ( ^OD to the Fourier transform of this equation, and substitute 



(0) and (ID to obtain 
from which it follows that 

T^Cp,m-p = — {l2[rn~p) ~ 7p ~ In-p) Cp,m-p- (55) 

We now invoke the standard solution of ([55|), applicable at the uv fixed point [12], to 
obtain the scaling relation 

Cp,m-p{r) ~ r-/3--, (56) 

where 

^np = 72{m-p) - 7p - 7^-p- (57) 
Upon substituting ( |56D in (|52D, it is immediately evident that the scaling exponent of Sn{r) 



is given by 



71 

o - Tn, (58) 



where 

Tn = max Tnp, for n = 2m > 2. (59) 

p 

Once 7* has been evaluated from (^), at the fixed point, which we do in Section |V], it is a 
simple matter to evaluate r„, as we show in Section 

Odd orders with n = 2m + 1 > 3 may be treated similarly with minor adjustments 
to allow for the fact that the expansions involve odd powers. In this case, however, it is 
immediately evident that the dominant scaling must arise from the Wilson coefficient of the 
unit operator corresponding to p = m, because averaging wipes out other terms by virtue 
of the fact that (O2S+1) = 0. Hence, we obtain 

rn = -{l*m + l*m+i) ^oi n = 2m + 1 > 3. (60) 



Again, the justification of the relevant expansions is given in Section |IX 
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IV. THE LINEAR RESPONSE 



In order to evaluate r^, we have to establish that an uv fixed point exists, which entails 
showing that the RG /5 function (^Tj) possesses a zero 

Pig*) = 0, (61) 

at which 

dp/dg < 0. (62) 

To do this we must first determine the dependence of the renormalisation constants and 



Zf) on the renormalisation scale /i. We will then verify that (|37D holds and use this fact to 
calculate g^, from Z^. 

Consider Z,^. According to the general theory of renormalisation [14], we have an expan- 
sion of the form 



Z^ = l + ga,^ log (9 +9'[^ W + log } + - " 



(63) 



Here k is the wavenumber cut-off which provides the intermediate regulation of the divergent 
integrals. This is an ir wavenumber of the order of L^^, where L is the typical length scale 
of the large scale flow. Divergences arise in the limit k ^ 0, corresponding to the inertial 
range limit r/L 0. The constants ai^, and will be calculated by eliminating the 
logarithmic divergences, at 1 and 2-loop orders respectively, from the IPI Green's function 
Ta/3{k, /), which is the inverse of the Fourier transform Gap{k, 1) of the linear response function 

Go.p{x.x) = (^^\. (64) 



,5//3(£') 

Tafii and the other IPI functions that we shall require, are generated from the functional 
which is obtained in the usual way by performing a Legendre transformation on Wc = 
logH^, with respect to the sources of the elementary flelds, J and J, while holding the 
composite operator sources tg flxed [7,13]. The new source flelds for K are therefore given 
by 
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and 



,4 ""c 



Dk. 



i63{-k) 

with X itself given in terms of its source fields by 

K{u,u,ts) = -Wc + i j |j(-A;) ■u(A;) + J(-A;) 
It follows, therefore, that 

T^f,{kJ) = (27r)8-— ^ ^. 

i6ua{k)6uf3{l) 

Introduction of the reduced forms 

G^p{k, I) = {2n)^sCk + l)Pap{k)GCk), (65) 

and 

T^fsikJ) = i27i)^5Ck + i)Pa(,{k)TCk), 
then leads to the standard relation 

r(fc) = GCk)-\ (66) 

We can now use ( [3l| ) to show that the connection between the bare and renormalised forms 
is 

Tr = Z^Tb, (67) 

which demonstrates the suitability of r(A;) as a basis for determining Z^,. 

In carrying out the renormalisation of r(A;) to obtain the coefficients in (|63|), we choose 
the normalisation point to be A; = m, where 

171 = (m, = 0). 
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Here m is a fixed vector of magnitude 



|m| = fi, 

the direction of which need not be specified, because the geometrical factor is contained in 
Papim) which cancels off. The expansion is used in conjunction with a normalisation 
condition that sets r(m) equal to its tree level value. Thus, from (|l^) and (|66D , we have 

r(m) = ro(m) = Go{m)-' = r(/x)-\ (68) 

and so the 1-loop term satisfies the normalisation condition 

ri(m) = 0. (69) 

The Feynman diagram giving the 1-loop term of Taisirh) is shown in Fig.3(z). The 
standard rules apply to such diagrams with the following assignments, which are shown in 
Fig.l: 

1. External lines represent functional differentiation with respect to u{k) when continu- 
ous, and u(A;), when dotted. The diagram is divided by a factor of i for each differen- 
tiation with respect to u. 

2. A continuous line linking two vertices denotes the reduced velocity correlation function 
defined through 



vMvfsil)) = {2Txf6{k + l)Q^p{k), 
and given by 

Q^p{k) = D{k)\G{k)\^ P^^ (k). 

For ease of notation, we omit zero-order labels in writing down mathematical expres- 
sions for the diagrams. 
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3. A half dotted/half continuous line connecting two vertices represents i times the zero- 
order response function 



4. The NS vertex with one dotted and two continuous lines represents Pai3j{k.), the ar- 
gument of which is associated with the dotted leg, with k directed away from the 
node. 

Returning now to the 1-loop diagram for Ta/3{m), we note that it has a symmetry factor 
of 1. Hence, it yields a contribution to ri(m) given by 

Pa/3(m)r'i(m) = J DpPajs{ra.)Pxuf3{ra.-p)Gjx{m-p)Qsu{p)- 

We can extract the logarithmic divergence from this integral by expanding its integrand in 
powers of p/m. This is possible because the divergence emanates from the region p k, 
while K <^ /i. A simple calculation leads to 

r'i(m) = Igri^^r'loie), (70) 

where 



oo 



W) = [ ^7^, (71) 

£ 

in which the lower limit of integration is 

r(«) ■ 

Extracting the logarithmic singularity from this integral gives 

r;(m) = -r(/.)-i5log(^). (72) 

To this we have to add the term arising from the counterterm vertex shown in Fig.2(z). This 
contributes the term P^p (m) r'l{m) where 
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r'/(m) = AZ^t{h) ^ = aiuT (/i) ^ glog (J^ 



(73) 



But, from the normalisation condition (169|), we have 



T[{m) + T'l{m) = 0, 
which, upon substituting ( |72D and (j7^), yields 

aij, = 1. 

We next carry out the analogous calculation for Zd and show that its corresponding 
coefficient Qid also equals 1, thereby verifying that the condition (^) is satisfied at 1-loop 
order. Here the relevant IPI function is the correlation function given by 

which is readily shown to be related to the velocity correlation function Qap{k, I) by [13] 



Tiaf3{k,l) = / T^x{k,p)T(^f,{l,q)Qx^{p,q)dpdq. 



Substituting the reduced forms 



n,^(fc,/) = {2nY6{k + i)Paf3{k)U{k), 



and 



we get 



n(^) = r(A;) QCk) 



From this result and (|31| ) and (pq), we find that the bare and renormalised forms of 11 are 
related by 



= ZdUb, 
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which confirms that n(A;) is the appropriate IPI function to use for calculating Zd. 

The normalisation condition is again chosen to be consistent with the tree level approx- 
imation. That is, we set 

n(m) = no(m) =D{i£), 

so that the 1-loop term satisfies the normalisation condition 

ni(m) = 0. (74) 

The 1-loop Feynman diagram for Jla^irn) is shown in Fig.3(M). It has a symmetry factor 
of 1/2, and makes a contribution to ni(rfi) which is given by 

P„/3(m)ni(m) = -Dp-Pa75(m)P^Ai.(m)Q^A(]3)Q5,.(m - p). 

In extracting the logarithmic singularity from this integral, we must take into account the 
fact that the symmetry of the integrand results in singularities of equal strength at both p ~ 
and |p — m| ~ /t, the effect of which compensates for the symmetry factor. Consequently, 
we get 

n;(m) = -Z}(/i)^7log(^). (75) 
The Zy) counterterm, which is shown in Fig.2(M), contributes a term to J\\{rn) given by 

n;'(m) = ^ZoD (/i) = ai^D (/i) g log . (76) 
But, from the normalisation condition (fflD, we have 

n;(m) + n;'(m) = o, 

and substitution of ( ffSD and (|76D leads to 

a\D = 1- 

We shall take the equality of the 1-loop coefficients of Z^, and Zd as establishing that 
(P7|) holds. This allows us to calculate the uv fixed point from the linear response function 
alone as follows. We use the standard result [14] 
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m = -9f^§j^iogZ,(^l + g^^\ogZ.)j , (77) 

where Zg is the renormahsation constant associated with the couphng constant. From the 
definition Zg = qq/q and (|6D,(|2|),(|3g) and {^), we get 



Inserting this result in ( [TTD and substituting the expansion (|63D, leads to 

^[g] = 2g\l + a2ug). 

This yields an uv fixed point 

g* = -—, (78) 

which satisfies (|6lD and (^) provided that < 0. It remains,then, to calculate a2u- 

The constant is obtained from the 2-loop term of Taiji^k), namely Pa/3{m)T2{k). At 
the normalisation point it must satisfy the condition 

r2(m) = 0, (79) 

by virtue of (^91). Only two Feynman diagrams yield logarithmic divergences. They are 
shown in Figs.3{iii) and (zf).They contribute the terms 

Pa/3(m)r2(m) = - j DpDqQseip)Qxpiq) 

xP«^5(m)P^^,(m - p)P^xu{-m. - p)Pg^^(m - p - q), 



and 



P,/3(m)r^'(m) = - J DpDqQs.{p)QxM (80) 
xG^^{m - p)Ggy{rh -p- q)G^r{m - q) 
xPa^s{m)P^^p{m - q)P«,Ai/(m - p)Pg^^(m - p - q). 
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To extract the logarithmic singularities from these integrals, we expand their integrands 
in powers of both p/m and q/m. We do this in two steps. First, we integrate over frequencies 
and solid angles to get 

V',{m) = -^/r(/x)-i/i (e) , (81) 



and 



where 



n{m) = -^-g'r{f,)-'h{e), (82) 



h {e) = [ [ , (83) 

J J x^y^(l + x) {1 + X + y) 

e e 



and 



f f dxdy 
^^^^ J J x^yHl + x){l + y){l + x + yy ^^'^^ 

e e 

Secondly, we expand these double integrals for small e to obtain 

r^(m) = 8^7M/^)-Mog (^) , (85) 



r'2W = ?^r(/i)-Vlogf^). (86) 



and 

2 '"^^ 

To these two contributions to r2(m), we must add the counterterm, which, by analogy 
with ([73|) , takes the form 

n"(m)=a2.A(/i)-'log(^). (87) 
Thus, the normalisation condition (^) becomes 

r'2(m)+r'2'(m) + r^"(m) = 0, 

which, by (^)- (|87|) , yields 
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37 

a2. - -y. 

Therefore, from ([78|), we obtain the fixed point couphng constant 

9* = ^, (88) 

which verifies that the residual couphng can be treated as weak. 

FinaUy, we explain why the 2-loop topologies, which have been discarded in calculating F, 
do not contribute to 02;/. As we shall explain further in Section |Vl| , divergences arise in these 



diagrams when it is possible for one or more soft wavevectors (ie values of p and/or g <^ m ) to 
flow through a correlator. However, if this entails the flow of some or all of these wavevectors 
through the active (ie dotted) leg of the NS vertex, then the logarithmic divergence will be 
suppressed by the extra powers of p and/or q. In the case of the 2-loop diagrams which 
we have just calculated, the external hard wavevector rh flows through the active legs of all 
vertices, so no suppression occurs. However, in the case of the remaining topologies at least 
one soft wavevector flowing through a correlator must also flow through the active leg of a NS 
vertex. In the case of the four remaining two loop topologies containing vertex corrections, 
the logarithmic divergence is suppressed individually for each diagram, after integration over 
the solid angles. In the case of the three remaining 2-loop diagrams containing insertions of 
the 1-loop diagrams {i) and {iv) of Fig. 3, suppression results after integrating over the solid 
angles and summing over the diagrams, the overall cancellation being related to the fact that 
the coefficients aiy and aio associated with the two types of insertion are equal. Likewise 
the four 1-loop diagrams containing the counterterm vertices yield no net contribution to 
a2jy. The treatment of the power and power x logarithmic divergences arising in integrals 
like (|7TD,(p3D and (Bl) is given in Section |VI4 For the moment we discard them because 



they are not directly relevant to the actual calculation of the scaling exponents for reasons 
already given. 
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V. THE NONLINEAR RESPONSE 



Having established that an uv fixed point exists, we can proceed with the calculation 
of the anomalous dimension 7^ of the general operator Os{x), which is required for the 
evaluation of the anomaly r„.To do this in the simplest possible way, we must identify a IPI 
response function which can be renormalised by means of Zg. Elimination of the logarithmic 
divergences from such a function will then enable us to determine the constants in the 
expansion 

z. = 1 + a A" log (9 +£{\ (4") 'log' © + 4" log (9 } + ■ ■ ■ ■ (89) 

so that we can calculate 7^ using (^3|). 

Consider first the case s = 2. Obviously, the required function must involve 02(5;), which 
is the composite operator associated with the longitudinal turbulence energy. In addition, 
it must involve the dynamic response operator ([Tl|) in order to relate the anomaly T2 to 
the dynamics of the turbulence. This suggests that we should consider how the turbulence 
energy responds on average to a change in the forcing. Clearly, we can characterise the 
response of the turbulence energy at a point x to a change in the forcing at two points x' 
and x" by means of the nonlinear Green's function 

But the complexity of this object is such that its logarithmic divergences cannot be summed 

using the renormalisation group in terms of the Z2 and Z,y counterterms alone. On the 

(2) 

other hand, its average G^i^{x) taken over the forcing separation x' — x", which gives a mean 
response to forcing at the centroid of the excitation points, can be, as we shall show shortly. 
Hence, its corresponding IPI function provides a direct means of obtaining the expansion 
(^) and so it provides an adequate basis for the calculation of 72. 

This IPI function is obtained as follows. We start with the Fourier transform of (pO[), 
the reduced form of which is given by 
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(91) 



where 



(92) 



Its corresponding IPI response function follows from 

5^K 



ei'j(fc,/,p) = (27r) 



12 



with a reduced form given by 



A standard calculation shows that it is related to G]^g by 



^alik, l,P) = -j ^xa{q, k) r^^(g', /) G'^^iq, q',p) dqdq' 
from which, on making use of (|9TD-(p^, we obtain 

e(2)(A;,/) = -r(A;)r(/)G2(A;,/). 



(93) 



(94) 



(95) 



(96) 



(2) 

Next, we average G)^^ over the forcing separation to get 



Gap{x) = 2 / Gi^j(A;, k) exp{2ik ■ x)Dk. 



This integral shows that the Fourier transform of G^p{x) depends only on the diagonal 
components of the reduced function (P^. It follows, therefore, from ( PB| ) and (0), that the 
IPI object which we need to consider, in order to determine Z2, is 

e(2)(A;,A;) = -V(kfG2{k,k). 



(97) 



Indeed, an application of (pi]), together with (|67D, shows that its bare and renormalised 
forms are connected by 



Qf{k,k)=Z2Q'B\k,k) 



(2), 



30 



In this way, as we have indicated, we arrive at a function which can be renormahsed using 
the Z2 counterterm alone. 

The normahsation condition for 0(^^(A;, A;) is again apphed at the point k = m, and 
chosen to be consistent with the tree level approximation, which gives 

e(^)(m,m) = ef,^^(m,m) = -1, (98) 

so that the 1 and 2-loop terms satisfy the normalisation conditions 

eS^^ (m, m) = ^ (m, m) = 0. (99) 

The diagrams giving 

e(2) 

(m, m) to 2-loop order are shown in Fig. 4. Their new feature is 
the appearance of the heavy dot vertex. This represents the O2 composite operator vertex, 
which is shown in Fig.l(zf) for the general case of Og.We can understand how these diagrams 
arise from the loop expansion of K by using the general procedure described in [15]. This 
depends on the fact that (0) is a special case of the Ath order correlation function of 
elementary fields defined by 

in which the arguments x and z coalesce. Hence, their Fourier transforms are related. In 
particular, the connection between their respective IPI functions is 

eJJ](A;,/,m) = ^apxf.ikj,m - q,q)Gxiim - q)G^i{q) Dq, 

where $ is the IPI form corresponding to B^^\ which is generated by 



^af3^s{k,l,p,q) = (2 



,16 



Sua{k)6uf3{l)i6uj{p)i6us{q) 
This implies that the diagrams for Q^'^\m,m) are constructed from the diagrams for $ by 
tying the two dotted external legs of the latter to form the O2 vertex. 

The 1-loop diagram for Q^'^\m, m) shown in Fig.4(zi') is constructed from the tree level 
diagram for $, which is shown opposite to it in Fig.4(z). Similarly, the two 2-loop diagrams 
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for Q^'^\rh,m) , shown in Figs.4(v) and (vi), are constructed from the 1-loop diagrams for 
$, again shown opposite to them in Figs.4(zz) and [iii). The other possible 2-loop diagrams 
for Q^'^\rh,m) , which arise from the two remaining 1-loop diagrams for $, are discarded 
because the logarithmic divergences disappear, after integration over the solid angles. In 
addition, diagrams which produce longitudinal terms obviously make no contribution to Z2 
and can also be discarded. 

The 1-loop diagram of FigA{iv) contributes to Q\ the term 

Poi(m)P/3i(m)ef^'(m,m) = J DpPx-ya{m-p)P^sf3iT^+p)Gxi{-m+p)Gui{-m-p)Q^s{p), 

which yields a logarithmic divergence 

ef^'im, m) = -g log (^^) . (100) 



The contractions implied in ( pS]) again permit us to discard the longitudinal part of the 
above integral. The counterterm vertex shown in Fig.2(iii) adds a contribution 

-P„i(m)P^i(m)ef)"(m,m) = -P,i (m) P^i (m) AZ^, 

SO that by (|^) its contribution to Gj^ is 

ef ^"(m, rh) = -af^g log . (101) 

But, from the normalisation condition (PU]), we have 

Qf\rh,rh) = Qf'^\rh,rh) + Q^p'\m,rh) = 0, 

which, upon substituting (|100|) and (|101|) , gives 



,(2) 



At 2-loop order the diagrams in Figs.4(t') and (vi) contribute the terms 

(m,m) = -g Jg, 
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(102) 



and 



^(2)/// . . N y 2 7- 



where 




2 + x + y 



x2y2(i + a;)(2 + y)(l + X + ?/) 



dxdy, 



and 




(2 + x)(2 + X + y)(l + 3?/ + - (1 + x)y{2 + y){3 + x + y) 



xV(l + x)(2 + + ?/)2(2 + + x + 
The latter yield logarithmic divergences 

ef'(A.A) = 59^(i-ilog2) log(£ 



dxdy. 



(103) 



and 



Gr(m,m) = ^/ (^ + ^1°S2) log (^) . 



(104) 



To these we must add the 2- loop counterterm corresponding to ( |101D , namely 

^(2)W/ ^ ^ N (2) 2 1 



But the normalisation condition ( p9D gives 



©f^'(m, m) + Oi^''"(m,m) + Q^^'''\m,m) = 0, 



(2)W/ 



(105) 



which, after substituting ( |103| )- (|105| ), yields 



a?^ = 7.0. 



(106) 



The foregoing can be generalised to arbitrary s. In place of (|93D, we now consider the 
general IPI response function 



©SU(^i,-,^s,p) = (2vr)^(^+^) 



5ua^{ki)...5ua,{ks)5ts{-p) 
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with a reduced form defined by 

Then Zg can be found by eliminating the logarithmic divergences from the diagonal compo- 
nent Q^'^\m, . . . , rh) as above. The relevant diagrams are again those shown in Figs.4(zi')-(m), 
except that the heavy dot now symbolises the Os vertex of Fig.l(zw), so the s — 2 external 
legs of Os are not shown explicitly. Each diagram has a symmetry factor s(s — l)/2. As this 
is the only respect in which these diagrams differ from those just considered, we have the 
relation 

41 = i^^og. (107) 

However, this is an approximate result, because it is not valid for diagrams containing more 
than 2-loops. But, as we discuss further below, it suffices for the calculation of low order 
exponents. Thus, we have now calculated all the numerical constants that we require for 
the evaluation of C,n- 

VI. THE SCALING EXPONENTS 



For n = 2, we have, from 



2 

C2 = 3 + A2, 



where, from (H3D,(|48|) and CBgl). 



^2 = -gM? +oS'^9*)- (108) 

Substituting the numerical values calculated above, as given in (p8D,( |102]) and ( [L06| ), we get 

46 

A2 = ^ = 0.0336, 
372 

which yields 



34 



(2 = 0.70. 

For 77, = 3, we shall verify in Section |V111| that the known exact result 



C3 = l, 

holds. 

In the general case, for n > 3, we have from ( 



C =--T 



For even orders n = 2m, the anomaly is given by (|59|), 

r„ = maxr„p. (109) 

p 

But, from (||) , (§7p , , ([TO^) and (|TUS|), we have 



T"np = {p{p - 1) + {n-p){n-p - 1) - 2{m - p)[{2{m - p) - 1]} — . 
A simple calculation shows that the maximum value of this expression is attained by the 



two terms in the series (|5^) with (a) p = m and (b) p = m — 1; which gives for ( 109 ) 

Tn = m{m - 1)A2. (110) 
For odd orders, n = 2m + 1, the anomaly is given directly by (|60D, which yields 

where we have again used (^),(p^),( p.07| ) and (|108|) . 

The above results have been used to calculate C„ up to n = 10. The results are shown 
in Fig.5, together with the experimental data taken from [16-20]. It can be seen that the 
agreement is good up to about n = 7 and fair beyond, if we allow for the uncertainties in 
the experimental data which begin to arise. In particular, it may be noted that the key 
values (2 = 0.70 and Ce = 1-8 are in good agreement with experimental data, the respective 
data sets from [16-20] giving for (2 the values (0.71, 0.70, 0.71, 0.70, 0.71) and for the 
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values (1.78,1.8,1.8,1.71,1.71). The divergence of the experimental data at higher orders 
reflects the fact that the experimental determination of (n is not yet fully satisfactory for the 
reasons given in [20]. Hence, the good agreement between our calculations at higher values 
of n with the particular data sets from [16-18] must be treated with caution, particularly 
as the expression we have derived above is not applicable at large orders. This limitation 
stems from the fact that the mean nonlinear response function, being an average over the 
forcing conflguration, does not represent the effect of multiple correlations with sufficient 
accuracy at large n. In addition, the approximation ( |107] ), as we have noted, only holds 
up to 2-loop order. Indeed, it is evident from the foregoing that the overall approximation 
must fail when ng^, ~ 1. However, this occurs at roughly n = 20, which is well above the 
current limit of reliable experimental data. Equally, the divergence of our theoretical values 
at higher values of n from the other two data sets [19,20] could indicate that the accuracy 
of our low order approximation is already beginning to deteriorate at around n ~ 10. 



VII. ELIMINATION OF SWEEPING 

We now return to the question of the power and power x logarithmic divergences which, 
up to this point, we have simply discarded. The fact that power divergences arise when 
field-theoretic methods are applied to turbulence, using an Eulerian approach, was noticed 
originally in [21]. Their origin was subsequently identified as being due to the kinematic effect 
of the sweeping of small eddies by large eddies, having an almost uniform velocity [22,23]. 
The remedy was to change from an Eulerian to a Lagrangian description, but this greatly 
complicates the subsequent analysis [24]. However, it has been shown that the elimination 
of sweeping can be accomplished more simply by transforming to a frame moving with the 
local velocity of the large scale eddies at some chosen reference point, [25,26]. We shall 
show that a similar approach can be used to eliminate the power and power x logarithmic 
divergences within the present framework. In this way, we shall demonstrate that, although 
we have started out from an Eulerian formulation, we ultimately obtain quasi-Lagrangian 
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approximations for the renormahsed functions. 

The problem, therefore, is to find a sweeping interaction term, AL^, say, which can be 
used to ehminate the effect of sweeping convection. To this end, we introduce a uniform 
convection U into W and average over its probabihty distribution, which we assume to be a 
Gaussian distribution oc exp(— t/^/2f/Q).This adds to L an additional interaction term given 
by 

772 /• ... 
ALu = — Y / ^ ■ m v(m) ■ v(-?7i) v(/) ■ v( - /) DlDm, (111) 

which represents the effect of a random Galilean transformation of the velocity field. We 
have not distinguished between v before and after the transformation for consistency with 
the earlier expressions, such as (0), and bearing in mind that the transformation does not 
affect the statistical averages required for the structure functions. 

To represent diagrammatically the additional terms which arise in the loop expansion 
of W after the inclusion of the sweeping interaction term we need to introduce a new 4-leg 
'sweeping' vertex of the type shown in Fig.6(«). The two wavevectors / and m in ( |111| ) 



enter this vertex along its continuous legs and leave along the dotted legs. A pair of legs 
carrying a particular wavevector must also carry the same vector index to represent the 
scalar product. Free wavevectors in a diagram containing one or more of these sweeping 
vertices are identified, as previously, by overall wavenumber conservation, together with 
conservation at any NS vertex. Each sweeping vertex then contributes a factor Uq\ ■ m, 
where 1 and m are the two wavevectors which enter the vertex along its two continuous legs. 
In all other respects the diagrams are to be interpreted in accordance with the rules given 
in Section |^ 

Consider now the set of diagrams, containing only NS vertices, which are associated 
with a particular Green's function or velocity correlator, Q, say. Let C^s denote any such 
diagram contributing to Q. We shall show that it is possible to generate all power and 
power X logarithmic divergences of any C^s from a single sweeping interaction of the form 



1 1 1|) . Let Cu denote any diagram containing at least one sweeping vertex. If Cu contains 
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no NS vertices at all, then it will only generate power divergences. But if it contains 
at least one NS vertex, it will also generate power x logarithmic divergences. The following 
topological argument demonstrates that the power divergences of Cns can be put into 1-1 
correspondence with the Cu diagrams relating to Q. 

Each factor t{k,) (or, equivalently, e~^) in a power divergence of C^s arises because it 
is possible for a soft wavevector q to flow through a particular velocity correlator without 
flowing through the active legs of the two NS vertices which it connects, as already discussed 



in Section This situation can be represented diagrammatically by contracting the corre- 
lator into a 4-leg vertex formed by merging the two NS vertices which it links, whilst leaving 
the hard lines in tact. This can be demonstrated as follows. First, the new vertex must 
consist of two in-coming full lines which carry hard wavevectors, 1 and m (say), and two 
outgoing dotted lines along which they leave. This is because two full legs disappear from the 
merged NS vertices and wavevectors leave NS vertices along the dotted leg. Furthermore, 
after integrating over the directions of the soft wavevector q, the two merged NS vertices 
generate, through contraction of the projectors, a factor proportional to the scalar product 
of the in-coming hard lines, 1 ■ m, while the two legs of a pair carrying the same wavevector 
acquire the same vector index. The final integration over the wavenumbers then produces 
the constant 



oo 

If 3 

— / q^D{q)T{q)dq = -gr{K)u\ 



So, such a vertex must, in fact, be of the sweeping convection type (|111|) , with a coefficient 
given by 

= ^gri^y, (112) 

which relates the rms velocity of the sweeping eddies to the strength of the nonlinear inter- 
action. Note that Uq is scale dependent, as it depends on the renormalisation scale fi through 
g and z/, and, hence, it differs according to the fluctuation scale on which the RG focuses. 
In physical terms, this reflects the fact that the rms velocity of the sweeping eddies depends 
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on the scale selected. 

Clearly, if a subset of correlators of Cats, each of which carries a soft wavenumber, is 
contracted into such vertices, in a manner which allows hard wavevectors to flow through 
Cms, then the result is a diagram which is identical to one of the Cu diagrams. Moreover, 
it is clear that there are always exactly as many ways to contract the correlators in C^s 
there are different Cu diagrams and that their symmetry factors must match. This argument 
demonstrates, therefore, the important point that the power and power x logarithmic diver- 
gences generated by the NS vertex must arise on account of the background of kinematic 
sweeping effects. Moreover, we also see that, in order to eliminate them, it is only necessary 
to introduce a sweeping interaction term into W of opposite sign to the one from which they 
can be generated, which, according to (|1 1 1|) and (|112|) yields the sweeping interaction term 



ALs = f 1 • m v(m) ■ v( - m) v(/) ■ v( - /) DlDm. (113) 

At{k) J 

Thus, the sweeping vertex shown in Fig.6(z) is taken to represent the algebraic factor 

3 

Vertex 6(z) = — gT{k)h'^\ ■ m. 
2 



Having inserted ( |113| ) into W one is then left with only the pure logarithmic divergences 



generated by the NS vertex, which, as we have shown, can be summed using the RG. This 
justifies our procedure whereby power and power x logarithmic divergences are discarded 
when calculating anomalous exponents. 

We now illustrate the cancellation of power and powerx logarithmic divergences in con- 
crete terms by eliminating them to 2-loop order from r(A;). This will demonstrate how the 
various symmetry factors match up. Consider first the 1-loop diagram for ri(A;) arising from 
the NS vertex. From our previous result ([70D, we find that its power divergence is given, at 
the normalisation point, by 

ri(diagram3(z)) = ^(744- (114) 
Here the Feynman rules applied to the sweeping vertex yield the single diagram of Fig.6(zz), 
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as we anticipate from the fact that the NS vertices in Fig.3(z) can be merged in only one 
way. In this trivial calculation yields 

3 t(k) 
Fi (diagram 6[iij) = —-g- 



which cancels (|114|), as required. 



Explicit verification that there are no power or power x logarithmic divergences in r(A;) 
at 2-loop order is less trivial. Consider first diagram (iv) of Fig. 3. The power divergences 
arising from this diagram follow from (^2|) which gives 

F2(diagram 3{iv) = |^ ~^ e ~^ -loge > . (115) 

For this diagram the corresponding sweeping diagrams are diagrams of Fig. 7. This 

follows from the Feynman rules and can be checked from diagram (iv) of Fig. 3 by first 
contracting its correlators individually and then together. By applying the Feynman rules 
to diagram (z) of Fig. 7 we obtain, at the normalisation point, 

Papini)T2 (diagram 7{i)) = ^gT{K)u^ j p • (k - p)PAai.(ni)P^p^(m - p) 



We can evaluate this integral using the method described in Section |^ This gives 



Fa (diagram 7(z)) = 7^^- / 
4r(/i) e J 



4r(/i) e J x^{x + 1) 

e 

9 g^ ( 1 2 , 1 , 

'- + -log£ + -^ 



4r(yu) [e^ e s 
Diagram (ii) of Fig. 7 yields the same value 

r2(diagram 7{ii)) = r2(diagram 7(z)). 

Finally, evaluation of the diagram [iii) of Fig 7 is trivial and yields 

9 g^ 1 

Fa (diagram 7 {in)) = 
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Evidently, the sum of these three diagrams cancels (|115|) exactly. 

Similarly, we can show that the sweeping vertex eliminates the power divergences arising 



from the second 2-loop diagram, shown in Fig.3(zM). From ([81|), these are given by 



9 ( 1 51 4 1 
r2(diagram = --—rr \ ^ + 7;- + - ^oge \ . (116) 

4r(yu) {e^ 2e e J 

In this case, the corresponding diagrams generated by the sweeping vertex are diagrams 

{iv)-(vi) of Fig. 7 which contribute the terms 

r2(diagram 7{ivJ) = |^ + ~^°S^J , 

r2(diagram 7iv)) = 

and 

r2(diagram 7{m)) = + ' l^g^ + 



Again, their sum exactly cancels (|116|) . We have thereby verified to 2-loop order that the 



sweeping interaction eliminates power divergences from the linear response function. 



VIII. THE KOLMOGOROV APPROXIMATION 

The fact that it has been possible to calculate the anomalies successfully by means of 
perturbation theory stems, in part, from the incorporation of the Kolmogorov theory into 
the zero order approximation. As we have seen, this has been done by replacing the actual 
viscous quadratic form in W, arising from the NS equations, by a modified quadratic form, 
characterised by an effective random stirring force spectrum D{k) and the effective timescale 
r(A;).We now demonstrate that these two functions can be deduced self-consistently as part 
of the calculation and confirm that that they do have the inertial range forms given in (p^ 
and (|2|). 

To determine these functions, we need two conditions. As in [10], one condition is 
supplied by evaluating the energy equation to 1-loop order, which gives the convergent DIA 
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form, corresponding to the so-called line renormalisation [24]. In the inertial range, it reduces 
to the condition that the energy flux across wavenumbers IlE{k) is independent of k and 
equal to the mean dissipation rate e: 

IlE{k) = e. 

Thus, evaluation of IlE{k) to 1-loop order gives the well-known result [24] 

oo 

UE{k) = j T{p)dp, 

k 

where 



T{p) = 871^ J J dqdr 



p^qr 



r(p)-i + r(g)-i + r(r)-i 
X {6(p,g,r)Q(r)(g(g) - Q{p)) + h{p,r,q)Q{q){Q{r) - Q{p))} . (117) 

Here A indicates integration over the region of the p, q plane in which p, g, r can form a 
triangle and 

ijp' + q^ — r^)^ — {p^ — g^)^ 



g, r) 



The second condition must be deduced from the linear response function. This is where 
difficulties have arisen with this approach in the past, when using an Eulerian framework, 
because of the ir divergences arising from sweeping. On the other hand, it is known that no 
divergence problems arise from sweeping convection in the case of the energy equation [24]. 
However, we have just shown how these power divergences can be systematically removed 
from the response function (and, indeed, all such functions) by means of a random Galilean 
transformation of the velocity field. This leaves the logarithmic divergences which, as we 
have seen, are to be eliminated from T{k) using the Zj^ counterterm. Recall that to fix the 
finite part of r(A;), after this renormalisation, we imposed the normalisation condition (68) 



which specifies that its tree level term should be exact at the normalisation scale /x.Thus, 



after eliminating sweeping convection, as described in Section VII, and using the 1-loop 
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normalisation condition (|68D to eliminate the logarithmic divergences, we obtain, at an 
arbitrary wavevector k (with a; = 0), the renormalised linear response function 

oo 

/i^r(/i)^ — /c^r(/c)^ f p'^T{pfD{p)dp 



i(k,U)-r(A:) + g^2^(^)^(^) j 





(r(A;) + r(p)) (r(/i) + r(p)) 
^"^T^iiY — k'^T{kY f p'^T{p)D{p)dp 



67r2 y (^(A;) + r(p)) (r(^) + r(p))' 



(118) 







It is precisely the condition that this expression should, indeed, yield a finite renormalised 
value which provides the required second relation, as we now explain. 

In the inertial range limit, we seek scaling solutions with r(fc) oc k~'^ and Q{k) oc k^, in 
which case D{k) = T{k)^^Q{k) oc k"-'^''. Now standard dimensional analysis shows that for 



( |117|) to hold in these circumstances, we must have a + 2b = —8, [24]. Furthermore, if this 
scaling solution were to produce a non-renormalisable divergence in the response function, 
it would arise in the second term of ( p.l8| ), since we can assume that a > 0. To prevent this 
from occurring, the coefficient of the integral must be zero, which requires 

r{k) ^ ///\2/3 
r(/i) \k) ' 

giving a = 2/3, and, hence, b = —11/3, so that a + b = —3. Thus, these relations do, in fact, 
yield the solution (123) and (P^), which we may conveniently re- write as 



r{k)-' = (3e^l^ei^ (119) 

and 

D{k) = ^e^^'k-^ (120) 
Therefore, the energy spectrum function 

E{k) = AnPQik) = 4nPD{k)T{k) 
takes the Kolmogorov inertial range form 
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and the integral in the third term of ( |118|) is, indeed, finite and yields 



T{k,0) = T{k)-'{l-g\og (-]}. 



For present purposes, explicit evaluation of the two constants is unnecessary, since they 
ultimately disappear from the calculation of the exponents, because they only occur through 
the coupling constant which, as we have seen, is eventually evaluated in terms of its fixed 
point value. 

Next, we comment briefiy on the effect of allowing for the perturbation terms (p3|) which 
give the difference between the modified quadratic form and the original viscous form. As 
in [10], we treat these terms as being of nominal order g. Their effect is, firstly, to re- 
introduce into r(A;) the viscous timescale t^, (k) which was replaced by TQ^k). Secondly, and 
more significantly, new divergences appear. However, it is not difficult to show that the 
divergent terms which are independent of h{k) and V sum exactly to the amount cancelled 
by the counterterms, as would be expected. In the inertial range limit V ^ 0, this leaves 
the term arising from which is given by 





Given that the actual stirring force spectrum function h{k) has remained arbitrary, subject 
only to the condition that it yields a finite input power given by 



it is clear that the above integral for AF must be finite. 

Thus, the role of these perturbation terms is not critical as regards calculating the 
anomalous exponents, provided that the the spectrum of the stirring forces is non-zero 
only at small k, as it should be. However, what we find is that, although forced at large 
scales, the above solution behaves in the inertial range as if the fiuid were stirred with a 
force spectral function oc A;~^. In this context, it is interesting to note that, in a study 
of the randomly forced NS equations by a stochastic force with zero mean and variance 



oo 




oo 
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oc [27], evidence of multiscaling of the structure functions has been found. In particular, 
the results obtained for the ratios Cn/C2 with the spectrum have been shown to agree 
with the values computed from the NS equations forced at large scales. This, of course, is 
exactly what one might expect from the above approximation. The present results are also 
consistent with the numerical calculations in [28], which suggest the scaling Ti^k) oc /c"^/^, 
as in (|119|) , for the Lagrangian micro timescale, as opposed to the scaling TE{k) oc k~^ for 



the Eulerian micro timescale, evidence for which has also been presented in [29]. As we 
have seen, the reason why the Lagrangian timescale applies in the present calculation is 
because we have eliminated sweeping by referring the velocity field to a frame moving with 
the local velocity of the large scale eddies which prevail at any chosen scale. This extracts 
the straining interactions, which shape the spectrum, from the background of convection, 
to yield quasi-Lagrangian approximations. 

In a sense, this derivation of the Kolmogorov quadratic form is analogous to a multiple 
timescale expansion in nonlinear wave theory, where part of the nonlinear behaviour is 
incorporated into the linear approximation, eg via a slowly changing wave amplitude, the 
variation of which is then determined from the nonlinear interaction by requiring the absence 
of secular terms in the higher order approximation. Here the requirement is similar in that 
it demands the absence of non-renormalisable terms in order to determine the nonlinear 
behaviour of the modified quadratic form. 

An integral part of the Kolmogorov theory is the exact result that in the inertial range 
limit 

S,{r) = -\eT, (121) 
5 

[5]. So we conclude this section by verifying that this result follows from the present treat- 
ment. 

Using standard symmetry relations, we can express 5*3 (r) in terms of the longitudinal 
component of the equal time triple velocity correlator 
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giving 

S,ir)=6Bni{r,0,0). (122) 
Now the general form of the Fourier transform of Bap-y must be 

and so F{k) can be expressed in terms of the transfer spectrum T{k) by 



F{k) = -T{k), 



while T{k) is given to 1-loop order by ( |117|) . Substituting these results in (|122|) gives 



Ss{r) = j ^fci - ^ exp(2fcir)Dk. 
This integral can be expanded in powers of r the lowest order term giving 



After integrating over the solid angle, we get 



CO 



S^{r) = -^j T{k)dk. 
ft 

This latter integral is, of course, the transport power IIe{h), which is a finite quantity at 
1-loop order and equal to the mean dissipation rate, as indicated in above, and, hence, we 



recover ( [121| ). 

The correlation function i?iii(x) also has an important role in the derivation of the OPEs 
required for the structure functions with higher odd orders, as we shall see shortly. 

IX. DERIVATION OF THE OPES 



We give finally the derivation of the dominant terms of the OPEs which we have used in 
Section |ITl| to obtain the structure function expansions. We deal first with the expansions 
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required for the higher order structure functions with orders n > 3. These can be obtained 
using the technique described in [30]. We defer discussion of the particular case n = 2 until 
last, because it requires a different approach for the reasons given in Section 

We begin by considering the OPE of the general product, defined in P9|), as it appears in 
the expansion (^) for S'„(r), taking first the case of even orders n = 2m, withp = 0, 1, . . . , m, 
namely 

n—p p 

A r ^ - ^+ ^- 

p\[n — py. 

where, as previously, f± = vi{x ± r /2,y, z,t), and we have used the definition (P). Let us 
consider the effect of inserting A„_p.p into a correlation function containing an arbitrary set 
of elementary fields Va^ix), . . . ,Vai{xi), as in (0). Then, following the approach of [30], we 
can derive the dominant terms which we have used in Section |IT| by considering how many 
of the f+ fields can be paired with a V- field to form products of lower order correlation 
functions. 

Consider the case p = m, ie 

{Va^ (Xl) ...Vat iXi)Am,mix,r)). 

Here each v+ can be paired with a V- to yield a product term 

{{V+V^)"") {Va^ (Xl) . . . (Xi)) , (123) 

which corresponds to the presence of a unit operator term in the OPE, [30]. If, instead, we 
only select m — 1 pairs of v^v_ products, we obtain a term of the type 




Now, in the limit as r — > 0, f+/2 behaves like an insertion of 02{x) into the correlation 
function of elementary fields [30]. Hence, this product tends to 

2 {{V+V.r-') {VaA^l) . ..VaAxi)02{x)) . (124) 



47 



But the averages of powers of simply yield non-stochastic functions of r, which we 

shall denote generically by Co(r), C2(r), . . . , as appropriate. Thus, from ( |123| ) and ( |124| ), we 
obtain, in the limit as r — > 0, 

Since the elementary fields are arbitrary, it follows that we have an OPE of the form 

Am,m(x, r) = Co{r)I + C2{r)02{x) + .... 

The point about expansions of this type is that the operators of increasing complexity do, 
indeed, produce subdominant terms in the expansion of Sn{r). Here, for example, the unit 
operator term, as we have shown, produces the dominant scaling with anomalous exponent 
given by (|110|) , whereas the quadratic term can be readily shown to give the smaller exponent 
Tn = [m(m — 1) — 1]A2, and, hence, is subdominant, while further terms in the expansion 
would produce even greater reductions. 

A similar argument applies when p = m — 1. In this case, however, we cannot pair every 
f+ with a f_. Therefore, the unit operator term cannot appear in the OPE for Am+i,m~i. If, 
however, we pair every f _ with a f + then the remaining v"^ pairs with the elementary fields 
and, in the limit as r — > 0, again appears as an 02{x) insertion. In this case, therefore, the 
OPE starts with 02{x) to give 

A^+i_^„i(x,r) = C2(r)02(x) + . . . . 

By continuing with this argument, we see that the dominant term of the OPE for the general 
case of An-p^p must take the form given in (|5l|). 

Consider next odd orders, n = 2m + 1. When p = m, we have a term of the form 

which, again, corresponds to the presence of a unit operator term, which is, thus, the 
dominant term of the OPE, giving 
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Am+i,m(^, r) = Co(r)/ + . . . . 
When p — m — 1, by pairing each i>_ with a we obtain a term of the form 

In the hmit as r — > 0, appears as an insertion of the cubic operator Os{x), so that here 
the OPE takes the form 

Am+2,m-i(£, r) = C3{r)03{x) + .... 
Continuing this process, we get for the next OPE 

Am+3,m-2(£, r) = C5{r)05{x) + . . . , 

and so on. But, in fact, the only term which contributes to Sn{r) for odd n is the unit oper- 
ator term of A^+i^^ because (02s+i(^)) = for any integer s, in the case of homogeneous 
isotropic turbulence. 

In the particular case of v+v_, we can establish the form of its OPE by using an expansion 
in the Fourier domain, in which the wavenumber q, corresponding to the separation r, tends 
to infinity, as described for instance, in [7,8]. To this end, we start by considering the general 
correlation function 

Ha!3\n{xi, X2 I x', x") = (Xi) Vp (£2) Vx {x') {x")) , 

for the case in which x' and x" tend to a common point x, well separated from xi and X2. 
For simplicity of presentation here, we have included only two arbitrary fields Va{xi) and 
V;3(a:2) -Denote its Fourier transform by 

Ha/3Xi.{p,p' I k, k') = (v^ ip) Vfs iP') vx{k)v^{k')) (125) 
= {2n)U{p + p' + k + k')H^px,.{P,P' I k,k')- 

Then, in terms of the reduced correlation function, we can write 
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p + p' ^ p + p' . 



/- p -\- p' 
DpDp'DqHai3x^{p,p' \ q —, -q 



X exp ^ip ■ [xi - ^ ^ ^ J +ip ■ (x2- ^ ^ ^ ) + iq ■ {x - x") \ . (126) 



When the arguments in ( p.25| ) coalesce to the common point x, we obtain the correlation 
function 

with Fourier transform 

Qaf^Xt.{P,P' I q) = {Va (P) ip') (vxVf,) (q)) 

= {'^T^f ^{P + P + q)Qap\^l{p^P I q)- 



Thus, corresponding to ( 126 ), we have 



Qap\f,{xi,X2 I x) = y DpDp'Qc,pxfji{p,p' I -P - p') exp {ip ■ {xi - x) + ip ■ {x2 - x)} . 

Let 'ifai3Xfi{p,p' I k,k') and 'E:a/3x^j.{p,p' \ q) be the IPI functions associated with the 
connected parts of Ha|3X^J.{p,p' I k,k') and Q{p,p' \ q). Denoting the connected part by 
superscript c, we have, as in Section |V|, 

Hi%,ip,P I k, k') = -Gaa^ (p) Gpp. ip') Gxx'Ck)G,,,Ck')^^,p,x','iP,P' I k, k'), (127) 

and 

Qalxi^iP^P I q) = -Gaa'{P)G(3p'{p')Ea'f3'X,,{P,P \ q) ■ (128) 

According to the standard procedure [7,8], the behaviour of the correlation function 
Hai3x^i{xi, X2 I x', x") as a function of x' — x", when x' and x" both tend to a common value 
X, can be deduced from the behaviour of '^apx^ip^P' I ^ — (p + ]3')/2, — g — (p + p')/2), in 
the limit as g ^ oo, which is, indeed, apparent from ( |126| ) and ( |127| ). Now, the diagrams 
which contribute to this IPI correlation function are diagram (z) of Fig.8, together with 
its permutation {p,a) ^ {P',P), and diagram (ii). However, it is easy to see from these 
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diagrams that, as g — oo, diagram (ii) yields a contribution which is smaller than that from 
diagram (z) by a factor Qap{q) ~ q^^^^^. So to derive the dominant term, we need to focus 
on diagram (z) and its permutation. The corresponding diagrams for Ea/BXfiipjp' \ —p — p') 
are diagram {in) of Fig.8 plus its permutation A ^ /i. 

Evaluation of these diagrams using the Feynman rules is straightforward and yields 

,Tr /I-- P + P' P + P\ fp / P + P\p / P + P\p / NO tn 
^apX,,{p,P I q ^) = / P^U^ ^)^Mr,y(-q —) Pa^^{p) Pf^usiP ) 

p' — p 

^Q^&{s)Qr^a{p - s)Qi^u{p' + s)Qpriq H h s)D? 

+ {p,a) + {p',p), 

and 



{p,p' \ -p-p') = - J Paay{p)Pl3u5ip')Q'y5is)Qaf,iP " s)Qxuip' + s)Ds 
+ (A ^ /i). 



Hence, for large q, we obtain from the last two equations the relation 
p ~\~ p' p -\- p' 

^af3X,^{p,p' I q -q —) = PxiM)Pt^Tn{^)Qpr{q)'^ap^^{p,p' I -P-P')- 



Combining this with ( |127|) and ( |128| ) yields the approximation 



rr(c) I . P^P' . P + ^ /| - /1on^ 

H'JxMp I ^ ^) = ^A^€'y(?)^AA^(^^'P I ), (129) 



where, to this order. 



Cx,Uq) = PxMP,rM \Gm'QA^)- (130) 



To obtain the required expansion for we must take the inverse Fourier transform 

of ( |129|) for the particular case A = yU = 1 with 



V T 

x' = {x + -,y,z,t) and x" = {x - -,y, z,t). 



The coefficient Ciigr,(x' — x") then depends only upon r and, according to ( |130| ), it must 
have the form 
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where F{q) is a function only of the wavenumber q. It is clear from this integral that Cn^^ 
must be diagonal in the indices ^,77, and have equal transverse components:Cii22 = Ciiss- 

We now define Q^^fji ) to be the connected correlation function formed from the 

elementary fields Va{xi) and f/3(x2), with the insertion of the longitudinal energy operator 
02(x), ie it is the particular case of (|36D with s = 2 and / = 2. Then 

<5Sll(^l>^2 I x) = 2Q^^^{xi,X2 I x). 

Similarly, we define Ql^f^ {xi,X2 \ x) to be the correlation function with Va{xi) and Vf3{x2), 
and the insertion of the transverse energy operator 

OP{x) = l{vl + vl). 

Thus, we have 

^1/322(^1' -^2 I x) ~^ Q^Jps^iXl, X2 I x) = 2Q^^p{xi,X2 \ x). 

Finally, we define longitudinal and transverse coefficients by writing 

C2(r) = 2Cnii(r), 

and 

C2(r) = 2Cn22(r) = 2Cn33(r). 

Using these definitions, and taking into account the diagonality of Cn^^^ enables us to express 
the inverse Fourier transform of ( |129| ), for the case A = = 1, as 

^a/]ii(*i'*2 I x,r) = C2(r)Qj^^^(xi,X2 I x) + C^(r)Q2^(xi,X2 | x), 
which, in the limit as r — 0, leads to 

{Vo,{Xi)Vp{x2)v+V_) = (va{.Xi)vp{x2) 
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E 



+ C2(r)02(:E) +C2(r)0f)(x) 



+ 



Since the fields Va{xi) and t'/3(x2) are arbitrary, we may conclude that 

E 

v+v. = -/ + C2(r)02(x) + .... 

Note that we have discarded the transverse operator because it is subdominant. This follows 
immediately from the analysis of Section 0. For example, in the case of when we 

/■ox 

calculate the corresponding value of the constant CL^ , clS defined in (g^), we get twice the 
value given in ( p.02| ) for the longitudinal operator 02(x), because, by isotropy, each of the 
two transverse components of 0^2^^ (x) contributes an amount equal to the value obtained for 
02(x) and, hence, the right hand side of (|108|) then yields an anomalous exponent of 2A2, 
indicating that 0^\x) makes a subdominant contribution to 5*2 (r). Thus, we have shown, 
to within the order g'^ of the calculation, that the dominant term of the OPE for v^v^ has 
the form given in (^). 

X. SUMMARY AND DISCUSSION 

The fact that it has been possible to demonstrate multiscaling and calculate anoma- 
lous exponents successfully from the generating functional by means of perturbation theory, 
notwithstanding the strong nonlinearity of the NS equations, is attributable to several fac- 
tors. These include: (1) the use of a modified quadratic form, which is derived self-consistently 
from the NS nonlinearity; (2) the incorporation in the generating functional of the composite 
operators which appear in the definition of the general structure function; (3) the application 
of OPEs to derive corrections to the Kolmogorov exponents in terms of the anomalous di- 
mensions of these operators; (4) the identification of a class of irreducible Green's functions 
containing insertions of these operators, which facilitate the calculation of their anomalous 
dimensions; (5) the elimination of sweeping convection effects using a random Galilean trans- 
formation of the velocity field; and, finally, (6) the deduction of the inertial range scaling 
using an uv fixed point of the RG to achieve the required small wavenumber limit. Let us 
now consider how each of these factors contributes to overcoming the obstacles encountered 
in previous applications of the RG. 

53 



The use of the modified quadratic form is an important element in the success of our cal- 
culation, because it provides an accurate initial approximation, which yields the Kolmogorov 
distribution in the inertial range limit. By contrast, in the early work which employed a field 
theoretic RG [31], and in subsequent developments of it [32-34], including equivalent formu- 
lations based on [35], reviewed recently in [36], the zero order approximation is based solely 
on on the linear terms of the NS equations, as in a conventional field theory calculation. 
Because this is a poor approximation for turbulence, it does not result in a genuine weak 
expansion parameter. For example, in the previous applications of RG techniques based on 
an expansion in the force spectrum exponent (ie the e-expansion) , in which the expansion 
about e = is extrapolated to e = 4, the value of the coupling constant is not small, at 
the ir fixed point which is used. Therefore, the accuracy of the expansion is uncontrolled. 
Indeed, according to [37], it may even be uncontrolled when e <^ 1, and there are problems in 
establishing its radius of convergence and the value of e at which long range driving becomes 
technically irrelevant [38]. 

However, our expansion is of a different nature. First, we do not use an e-expansion. 
Actually, there is no force power spectrum in our calculation as such. As we showed, the 
force spectrum h{k) remains in the calculation as an arbitrary function, subject only to the 
requirment that it yields a finite input power. What the modified quadratic form provides, 
however, is an apparent force power spectrum D{k), but its exponent is fixed by the solution 
( |12CI|) , and, thus, cannot be varied. Second, we do not use an ir fixed point, because we are 
interested in taking the short wavelength limit, for which purpose we require an uv fixed 
point. Together, these differences result in a genuinely small coupling constant g, which is 
about 1/20 at the fixed point, as shown in Section |^ Hence, our expansion is inherently 
more accurate than the e-expansion. In fact, given that our calculation is carried out to 
2- loop order, its errors are controlled at ~ 10^^. Another significant consequence of using 
the modified quadratic form is that no convergence problems are encountered in the uv 
region. This, together with the fact that we do not use an e-expansion or an ir fixed point, 
means that none of the ingredients which cause marginality by power counting in previous 
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applications of the RG [37], are present in our approach. 

On the other hand, there is a similar problem to be faced in the present calculation. Any 
fully renormalised theory of turbulence must contain an infinite number of renormalised 
functions because it must be equivalent to the hierarchy of equations for the cumulants.This 
equivalence has been demonstrated recently [39]. In fact, each cumulant will have a repre- 
sentation as a expansion in terms of irreducible renormalised functions. Thus, one has an 
infinite set of vertex functions to contend with. Now, when any one of these irreducible 
functions is calculated in perturbation theory using the modified quadratic form, the overall 
logarithmic divergence will remain, after sweeping divergences have been eliminated. So the 
problem in the present approach amounts to the resummation of these logarithms. However, 
we showed in Section 0that this difficulty could be overcome, in relation to multiscaling, by 
identifying the infinite sub-class of functions which yields the desired information relating 
to anomalous exponents while being, at the same time, amenable to resummation using the 
RG.The irreducible inserted nonlinear Green's functions defined in Section |V| satisfy both 
requirements. Being fully irreducible they give full n-point correlations. However, as we have 
seen, to render them tractable, it was expedient to obtain a mean response to forcing at the 
centroid of the excitation points. This averaging thus constitutes a closure approximation. 
Although this type of closure approximation permits considerable progress to be made with 
the calculation of the exponents, the averaging process limits its applicability to relatively 
low orders, n < 10, because the multiple correlations between the apparent forcing at differ- 
ent space-time points are not then approximated accurately enough at higher orders. Thus, 
a different approximation would be required to obtain the asymptotic scaling at large orders 
and it remains for future work to discover a suitable approach. 
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CAPTIONS TO FIGURES 

1. FIG.l: Components of the diagrams: (i) velocity correlator; (ii) linear response func- 
tion; (iii) Navier-Stokesvertex; (iv) composite operator Og vertex. 

2. FIG. 2: Counterterm vertices associated with the renormalization of the elementary 
fields and the compositeoperators. 



3. FIG. 3: The IPI Feynman diagrams for the linear response evaluated in Section IV 



4. FIG. 4: The IPI diagrams for the nonlinear response functions evaluated in Section 0. 

5. FIG. 5: Comparison of the theoretical expression for C,n (full line) with experimental 
data. 

6. FIG. 6: (i) The 'sweeping' vertex; (ii) the 1-loop 'sweeping' diagram for the linear 
response function evaluated in Section |V14 

7. FIG. 7: The 2-loop 'sweeping' diagrams for the linear response function evaluated in 
Section IVll. 



8. FIG. 8: The 1-loop diagrams for the correlation functions evaluated in Section [IX 
connection with the OPEs. 
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